Limits on Neutrino Emission from Gamma-Ray Bursts with the 40 String IceCube 

Detector 



R. Abbasi, 28 Y. Abdou, 22 T. Abu-Zayyad, 33 J. Adams, 16 J. A. Aguilar, 28 M. Ahlers, 32 K. Andeen, 28 

J. Auffenberg, 38 X. Bai, 31 M. Baker, 28 S. W. Barwick, 24 R. Bay, 7 J. L. Bazo Alba, 39 K. Beattie, 8 
J. J. Beatty, 18 ' 19 S. Bechet, 13 J. K. Becker, 10 K.-H. Becker, 38 M. L. Benabderrahmane, 39 S. BenZvi, 28 
J. Berdermann, 39 P. Berghaus, 28 D. Berley, 17 E. Bernardini, 39 D. Bertrand, 13 D. Z. Besson, 26 D. Bindig, 38 
M. Bissok, 1 E. Blaufuss, 17 J. Blumenthal, 1 D. J. Boersma, 1 C. Bohm, 34 D. Bose, 14 S. Boser, 11 O. Botner, 37 
J. Braun, 28 A. M. Brown, 16 S. Buitink, 8 M. Carson, 22 D. Chirkin, 28 B. Christy, 17 J. Clem, 31 F. Clevermann, 20 
S. Cohen, 25 C. Colnard, 23 D. F. Cowen, 36 * 35 M. V. D'Agostino, 7 M. Danninger, 34 J. Daughhetee, 5 J. C. Davis, 18 

C. De Clercq, 14 L. Demirors, 25 O. Depaepe, 14 F. Descamps, 22 P. Desiati, 28 G. de Vries-Uiterweerd, 22 T. DeYoung, 36 
O ; J. C. Diaz-Velez, 28 M. Dierckxsens, 13 J. Dreyer, 10 J. P. Dumm, 28 R. Ehrlich, 17 J. Eisch, 28 R. W. Ellsworth, 17 

O. Engdegard, 37 S. Euler, 1 P. A. Evenson, 31 O. Fadiran, 4 A. R. Fazely, 6 A. Fedynitch, 10 T. Feusels, 22 
K. Filimonov, 7 C. Finley, 34 T. Fischer- Wasels, 38 M. M. Foerster, 36 B. D. Fox, 36 A. Franckowiak, 11 R. Franke, 39 
T. K. Gaisser, 31 J. Gallagher, 27 M. Geisler, 1 L. Gerhardt, 8 - 7 L. Gladstone, 28 T. Gliisenkamp, 1 A. Goldschmidt, 8 
J. A. Goodman, 17 D. Grant, 21 T. Griesel, 29 A. Grofi, 16 - 23 S. Grullon, 28 M. Gurtner, 38 C. Ha, 36 A. Hallgren, 37 
F. Halzen, 28 K. Han, 16 K. Hanson, 13 ' 28 D. Heinen, 1 K. Helbing, 38 P. Herquet, 30 S. Hickford, 16 G. C. Hill, 28 

K. D. Hoffman, 17 A. Homeier, 11 K. Hoshina, 28 D. Hubert, 14 W. Huelsnitz, 17 J.-P Hiilfi, 1 P. O. Hulth, 34 
K. Hultqvist, 34 S. Hussain, 31 A. Ishihara, 15 J. Jacobsen, 28 G. S. Japaridze, 4 H. Johansson, 34 J. M. Joseph, 8 
K.-H. Kampert, 38 A. Kappes, 9 T. Karg, 38 A. Karle, 28 J. L. Kelley, 28 N. Kemming, 9 P. Kenny, 26 J. Kiryluk, 8 * 7 
F. Kislat, 39 S. R. Klein, 8 - 7 J.-H. Kohne, 20 G. Kohnen, 30 H. Kolanoski, 9 L. Kopke, 29 S. Kopper, 38 D. J. Koskinen, 36 
M. Kowalski, 11 T. Kowarik, 29 M. Krasberg, 28 T. Krings, 1 G. Kroll, 29 K. Kuehn, 18 T. Kuwabara, 31 M. Labare, 14 
S. Lafebre, 36 K. Laihem, 1 H. Landsman, 28 M. J. Larson, 36 R. Lauer, 39 R. Lehmann, 9 J. Liinemann, 29 
J. Madsen, 33 P. Majumdar, 39 A. Marotta, 13 R. Maruyama, 28 K. Mase, 15 H. S. Matis, 8 K. Meagher, 17 'B 
M. Merck, 28 P. Meszaros, 35 ' 36 T. Meures, 1 E. Middell, 39 N. Milke, 20 J. Miller, 37 T. Montaruli, 28 ^ R. Morse, 28 
S. M. Movit, 35 R. Nahnhauer, 39 J. W. Nam, 24 U. Naumann, 38 P. Niefien, 31 D. R. Nygren, 8 S. Odrowski, 23 
A. Olivas, 17 M. Olivo, 37 ' 10 A. O'Murchadha, 28 M. Ono, 15 S. Panknin, 11 L. Paul, 1 C. Perez de los Heros, 37 
J. Petrovic, 13 A. Piegsa, 29 D. Pieloth, 20 R. Porrata, 7 J. Posselt, 38 P. B. Price, 7 M. Prikockis, 36 G. T. Przybylski, 8 
K. Rawlins, 3 P. Redl, 17 E. Resconi, 23 W. Rhode, 20 M. Ribordy, 25 A. Rizzo, 14 J. P. Rodrigues, 28 P. Roth, 17 
F. Rothmaier, 29 C. Rott, 18 T. Ruhe, 20 D. Rutledge, 36 B. Ruzybayev, 31 D. Ryckbosch, 22 H.-G. Sander, 29 
M. Santander, 28 S. Sarkar, 32 K. Schatto, 29 T. Schmidt, 17 A. Schoenwald, 39 A. Schukraft, 1 A. Schultes, 38 
O. Schulz, 23 M. Schunck, 1 D. Seckel, 31 B. Semburg, 38 S. H. Sep, 34 Y. Sestayo, 23 S. Seunarine, 12 A. Silvestri, 24 
A. Slipak, 36 G. M. Spiczak, 33 C. Spiering, 39 M. Stamatikos, 18 'E] T. Stanev, 31 G. Stephens, 36 T. Stezelberger, 8 
R. G. Stokstad, 8 S. Stoyanov, 31 E. A. Strahler, 14 T. Straszheim, 17 G. W. Sullivan, 17 Q. Swillens, 13 H. Taavola, 37 
I. Taboada, 5 A. Tamburro, 33 O. Tarasova, 39 A. Tepe, 5 S. Ter-Antonyan, 6 S. Tilav, 31 P. A. Toale, 36 S. Toscano, 28 

D. Tosi, 39 D. Turcan, 17 N. van Eijndhoven, 14 J. Vandenbroucke, 7 A. Van Overloop, 22 J. van Santen, 28 M. Vehring, 1 
M. Vogc, 23 B. Voigt, 39 C. Walck, 34 T. Waldenmaier, 9 M. Wallraff, 1 M. Walter, 39 C. Weaver, 28 C. Wendt, 28 

S. Westerhoff, 28 N. Whitehorn, 28 'B K. Wiebe, 29 C. H. Wiebusch, 1 D. R. Williams, 2 R. Wischnewski, 39 
H. Wissing, 17 M. Wolf, 23 K. Woschnagg, 7 C. Xu, 31 X. W. Xu, 6 G. Yodh, 24 S. Yoshida, 15 and P. Zarzhitsky 2 

(IceCube Collaboration) 

1 III. Physikalisches Institut, RWTH Aachen University, D-52056 Aachen, Germany 
2 Dept. of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 354-87, USA 
3 Dept. of Physics and Astronomy, University of Alaska Anchorage, 
3211 Providence Dr., Anchorage, AK 99508, USA 
4 CTSPS, Clark-Atlanta University, Atlanta, GA 30314, USA 
5 School of Physics and Center for Relativistic Astrophysics, 
Georgia Institute of Technology, Atlanta, GA 30332, USA 
6 Dept. of Physics, Southern University, Baton Rouge, LA 70813, USA 
7 Dept. of Physics, University of California, Berkeley, CA 94720, USA 
8 Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 
9 Institut fur Physik, Humboldt- Universit'dt zu Berlin, D-12489 Berlin, Germany 
10 Fakultat fur Physik & Astronomie, Ruhr-Universitdt Bochum, D-44 780 Bochum, Germany 
11 Physikalisches Institut, Universitdt Bonn, Nussallee 12, D-53115 Bonn, Germany 



Oh 
6 



2 



12 Dept. of Physics, University of the West Indies, 
Cave Hill Campus, Bridgetown BB11000, Barbados 
13 Universite Libre de Bruxelles, Science Faculty CP230, B-1050 Brussels, Belgium 
14 Vrije Universiteit Brussel, Dienst ELEM, B-1050 Brussels, Belgium 
15 Dept. of Physics, Chiba University, Chiba 263-8522, Japan 
16 Dept. of Physics and Astronomy, University of Canterbury, Private Bag 4800, Christchurch, New Zealand 
17 Dept. of Physics, University of Maryland, College Park, MD 20742, USA 
18 Dept. of Physics and Center for Cosmology and Astro-Particle Physics, 

Ohio State University, Columbus, OH 43210, USA 
19 Dept. of Astronomy, Ohio State University, Columbus, OH 43210, USA 
20 Dept. of Physics, TU Dortmund University, D-44221 Dortmund, Germany 
21 Dept. of Physics, University of Alberta, Edmonton, Alberta, Canada T6G 2G7 
22 Dept. of Subatomic and Radiation Physics, University of Gent, B-9000 Gent, Belgium 
23 Max-Planck-Institut fur Kernphysik, D-69177 Heidelberg, Germany 
24 Dept. of Physics and Astronomy, University of California, Irvine, CA 92697, USA 
25 Laboratory for High Energy Physics, Ecole Polytechnique Federale, CH-1015 Lausanne, Switzerland 
26 Dept. of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA 
27 Dept. of Astronomy, University of Wisconsin, Madison, WI 53706, USA 
28 Dept. of Physics, University of Wisconsin, Madison, WI 53706, USA 
29 Institute of Physics, University of Mainz, Staudinger Weg 7, D-55099 Mainz, Germany 
30 Universite de Mons, 7000 Mons, Belgium 
31 Bartol Research Institute and Dept. of Physics and Astronomy, 
University of Delaware, Newark, DE 19716, USA 
32 Dept. of Physics, University of Oxford, 1 Keble Road, Oxford 0X1 3NP, UK 
33 Dept. of Physics, University of Wisconsin, River Falls, WI 54022, USA 
34 Oskar Klein Centre and Dept. of Physics, Stockholm University, SE- 10691 Stockholm, Sweden 
35 Dept. of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802, USA 
36 Dept. of Physics, Pennsylvania State University, University Park, PA 16802, USA 
37 Dept. of Physics and Astronomy, Uppsala University, Box 516, S-75120 Uppsala, Sweden 
38 Dept. of Physics, University of Wuppertal, D-42119 Wuppertal, Germany 
39 DESY, D-15735 Zeuthen, Germany 

IceCube has become the first neutrino telescope with a sensitivity below the TeV neutrino flux 
predicted from gamma-ray bursts if GRBs are responsible for the observed cosmic-ray flux above 
10 18 eV. Two separate analyses using the half-complete IceCube detector, one a dedicated search 
for neutrinos from ^-interactions in the prompt phase of the GRB fireball, and the other a generic 
search for any neutrino emission from these sources over a wide range of energies and emission times, 
produced no evidence for neutrino emission, excluding prevailing models at 90% confidence. 



Gamma-ray Bursts (GRBs) have long been proposed 
[H as one of the most plausible sources of the highest en- 
ergy cosmic rays, as the observed flux can be entirely 
explained if the primary engine of the bursts acceler- 
ates protons and electrons with comparable efficiencies. 
The electrons would produce the observed gamma-ray 
emission by synchrotron emission and, possibly, inverse 
Compton scattering, while the protons escape to form 
the high-energy cosmic rays observed at Earth. Waxman 
and Bahcall observed^ that, in this potentially 
detectable flux of high-energy neutrinos is produced by 
interactions when protons and photons coexist in the 
primary fireball. The detailed flux predictions are depen- 
dent on the fireball parameters; here we use the model by 
Guetta et al. 0] to compute these parameters from obser- 
vations by gamma-ray telescopes. Past searches with Ice- 
Cube and other neutrino telescopes have met with nega- 
tive results UHl], but have never before had sensitivities 
at the level of the expected flux. We search in this work 
for neutrinos in coincidence with 117 GRBs with half 
of the IceCube detector complete, and for the first time 



reach a sensitivity that would yield a positive result given 
expected fireball parameters, with a 4<r expected excess. 

IceCube is a TeV-scale neutrino telescope currently un- 
der construction at the South Pole which detects neu- 
trinos by measuring the Cerenkov light from secondary 
charged particles produced in neutrino-nucleon interac- 
tions. A total of 5160 Digital Optical Modules[7] contain- 
ing 10-inch photomultipliers and arranged in 86 strings 
frozen in the ice will make up the full detector; the results 
presented here were obtained using the first 40 of these 
strings. Although capable of detecting multiple flavors of 
neutrinos from the entire sky, for point sources the detec- 
tor is sensitive primarily to up-going muons produced in 
muon neutrino charged-current interactions. Searches in 
the muon channel benefit from good angular resolution 
(~ 0.7° for E v > 10 TeV) and from the long range of 
high energy muons (several km at TeV energies), which 
substantially increases the effective volume of the detec- 
tor. By using up-going tracks, the Earth is used to shield 
against the much larger flux of down-going muons from 
cosmic ray interactions in the atmosphere. Backgrounds 
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from cosmic ray-produced muons and atmospheric neu- 
trinos can be further reduced using the muon energy, as 
neutrinos from GRBs are expected to have higher ener- 
gies than from either atmospheric source. 

The origin of observed events in IceCube is determined 
by fitting a track to the hit pattern of the detected 
Cerenkov light using a maximum likelihood method [8]. 
The location of the maximum is used as the source of 
the associated neutrino (collinear with the muon), and 
the statistical uncertainty in the fit provides an estimate 
of the uncertainty on the reconstructed direction [9]. 

Due to the stochastic nature of muon energy-loss pro- 
cesses and the rarity of events fully contained within the 
detector, it is not possible to measure the energy of either 
the muon or the primary neutrino directly. It is, however, 
possible to measure the mean energy-loss rate of muons in 
the detector, which is correlated at high energies with the 
muon energy and with the original neutrino energyfioj]. 
The uncertainty of the muon energy using this method 
is on the order 0.3-0.4 in log 10 E. 

IceCube operated in a 40-string configuration from 
April 5, 2008 until May 20, 2009. During that time 
129 GRBs were reported in the northern hemisphere via 
the GRB Coordinates Network^ (GCN). We assembled 
a catalog using data from GCN notices and circulars, 
where the position of the burst was taken from the no- 
tice with the lowest reported positional error. For bursts 
which were localized only by the Fermi Gamma-ray Burst 
Monitor (GBM), the position was instead taken from the 
GBM Burst Catalog[12j. The start and stop times of the 
prompt gamma-ray phase, T sta rt and T stop respectively, 
were determined by taking the earliest and latest times 
any satellite reported detecting gamma-rays. The fluen.ce 
and spectral information were taken preferentially from 
Fermi GBM, Konus-Wind, Suzaku WAM, then Swift. 

Fermi GRBs for which no fluence was reported be- 
cause the burst was too weak were removed. GRB080521 , 
GRB081113 and GRB090515 occurred during detec- 
tor downtime and were removed from the catalog. 
GRB090422 and GRB090423 occurred during a prelimi- 
nary run with 59 strings in operation and will be analyzed 
later. The final catalog contained 117 bursts. 

Neutrino spectra were calculated [!, Q using data from 
the gamma-ray spectra of individual bursts, or average 
parameters if no spectral measurements were available. 
Definitions of parameters and equations used to calcu- 
late neutrino fluence are identical to Appendix A of 
Spectra were calculated as power laws with two breaks: 
a low energy break associated with the break in the pho- 
ton spectrum, and a high energy break from synchrotron 
losses of muons and pions (Fig. Q}. 

From the length of gamma emission and energy spec- 
trum, bursts are classified by GCN into two groups (long- 
soft and short-hard), which may have different underly- 
ing sources. If a burst was not explicitly identified as 
one class in a GCN notice, we used average values for a 



short-hard burst if 90% of the gamma emission was in 
less than 2 seconds and a long-soft burst otherwise. 
Parameters for average long-soft bursts are from 0| ■ For 
short-hard bursts, we used L™° — 10 51 erg/s, t yat = 0.001 
s, and for redshift (z) the average of all Swift short burst 
measurements. 
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FIG. 1. The neutrino spectra, including oscillations, of the 
five brightest GRBs are shown along with eight randomly se- 
lected bursts (thin lines). A single burst with Waxman 2003 
parameters fl3T] . assuming a cosmic ray energy density of 10 44 
erg Mpc -3 yr _ , is shown by a thin dashed line. The sum of 
all 117 individual bursts is shown as a thick solid line along 
with the Waxman 2003 p^] prediction in a thick dashed line. 

Two independent searches were conducted: one search- 
ing for neutrinos with the specific energy spectrum pre- 
dicted by Guetta et al. [3j during the period of maximum 
gamma emission, and the other searching generically for 
high-energy neutrinos within up to 24 hours of the ob- 
served bursts. 

The first of the two analyses, the model-dependent 
analysis, was designed specifically to find neutrinos pro- 
duced in interactions during the prompt phase of the 
GRB. Events observed in the detector were reduced by a 
series of cuts designed to select neutrino-like events, re- 
sulting in a data sample of primarily atmospheric neutri- 
nos, an irreducible background for this analysis. We then 
conducted an unbinned maximum likelihood search^] in 
which each event passing these cuts was assigned likeli- 
hoods of being a signal event (from a GRB) and of being 
a background event. Both the signal and background 
likelihoods for each event i were the product of three in- 
dependent probability density functions (PDFs) based on 
direction, arrival time, and muon energy. 

The spatial signal PDF was a 2-dimensional Gaussian: 

P sm = I expf- l fGRB - f *' 2 ^ fi) 

F [Xl) ~ 27 r(4 RB + af) ° XP I. 2(^ RB + af) ) W 

where |a?GRB — \ is the angle between the reconstructed 
neutrino direction and the best location of the gamma 
ray burst provided by GCN, and ctgrb and Cj are the 
localization uncertainty of the GRB and the muon recon- 
struction respectively. The spatial background PDF was 
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computed using a smoothed histogram of all off-source 
data in detector coordinates, accounting for zenith and 
azimuth asymmetry in the detector. 

The temporal signal PDFs were constant during the 
prompt phase of the Gamma-Ray Burst (between T sta rt 
and T s top), with Gaussian tails of width T st0 p — T s tart 
(constrained to at minimum 2 seconds and at maximum 
30). The background PDFs were constant in time. 

The signal energy PDF was computed from the recon- 
structed muon energy- loss (dE/dx) for neutrinos simu- 
lated with the average of the individual burst spectra 
(Fig. [T}, while the background energy PDF was com- 
puted from the dE/dx distribution of off-source data. 

From these likelihoods, we then computed the maxi- 
mally likely number of signal events. The resulting like- 
lihood ratio (the test statistic) was then compared to the 
distribution from scrambled background datasets to com- 
pute the significance of a result. 

As well as looking for neutrinos with properties mod- 
eled from measured burst parameters, we conducted an 
additional search (the model-independent analysis) using 
wider time search windows and looser event selection cri- 
teria, allowing observation of events with late or early 
arrival times or with unexpected energies due to unan- 
ticipated emission mechanisms. 

Starting at the interval from -10 to +10 seconds from 
the GRB trigger time, we expanded a search time window 
in one second increments in both directions out to ± 
one day, looking for a significant excess of neutrinos at 
each iteration. High correlation between adjacent time 
windows reduces the trials correction to the significance 
of any excess to only a few hundred. 

Event selection for the model-independent search was 
based entirely on rejecting misreconstructed downgoing 
atmospheric cosmic ray muons, which are the domi- 
nant background to this analysis, constituting more than 
99.9% of the final 161 million event sample. To avoid 
assuming a signal neutrino spectrum, no attempt was 
made to reject the small low-energy background from at- 
mospheric neutrinos. 

To ensure that no events were missed due to incorrect 
assumptions, this analysis was designed to maximize the 
number of signal neutrinos in the final analysis instead 
of the significance of an excess. Instead of being selected 
by hard cuts, events were weighted by their probabilities 
of being signal neutrinos Each probability was the 
product of the event's point-spread probability density 
function (Eq. [TJ and the probability that the event was 
a neutrino, determined by dividing smoothed histograms 
of detector data and neutrino simulation in several vari- 
ables related to reconstruction accuracy. These were then 
summed in each time window to form the expectation 
of the on-source signal neutrino density, which was then 
compared to the expected background value obtained by 
scrambling the observed data in time. 

Although the use of scrambled data for the background 



reduces many possible uncertainties, the use of simula- 
tion for signal introduces some systematic errors. The 
dominant sources of uncertainties in the final limits from 
both analyses are photon propagation in the ice, the 
quantum efficiency of the PMTs, and theoretical uncer- 
tainties in both the neutrino-nucleon cross section and 
cross-sections for muon energy-loss processes at high en- 
ergies. Depending on the analysis and time interval, the 
cumulative effect of these uncertainties amounts to 2-13% 
and has been included in the final limits using a Bayesian 



marginalization procedure [15 1. 

No events were observed in the model-dependent 
search with a signal to background likelihood ratio 
greater than one, with 2.99 signal events expected on a 
background of 0.097. The closest event to its associated 
GRB was 26° from GRB090301A. This sets a 90% upper 
limit of 82% of the expected flux in the region 37 — 2400 
TeV where 90% of the events were expected, including a 
systematic uncertainty of ~ 2% (Fig. [2]). 
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FIG. 2. 90% CL Neymanjlg] upper limit (including system- 
atics) set by model-dependent analysis in solid black with the 
expected Guetta et al. flux in dotted black. 22-string IceCube 
limit H is in dark gray and AMANDA in light gray. The 
Waxman 2003 flux[13| is shown for comparison in dotted light 
gray. Diffuse fluxes were obtained from fluences assuming a 
total of 667 uniformly distributed bursts per year. Fluences 
are aggregate for 117 bursts. 

In the model-independent search, no candidate events 
were observed in the interval ±2248 seconds with 4.2 ex- 
pected from the Guetta et al. calculation. The varia- 
tion of the upper limit (Fig. [3]) with At reflects statis- 
tical fluctuations in the background, as well as the pres- 
ence of individual events of varying quality. The three 
most significant of these occurred at —2249, —3594, and 
—6430 seconds respectively, and were low energy (~ 1 
TeV) neutrinos consistent with the atmospheric neutrino 
background. In addition to a constant 1|% uncertainty 
on the effective area (the ratio of fluence to the expected 
number of events), there is a systematic uncertainty in 
the limit on the number of expected events that increases 
with the size of the time window from 0-10% (included in 
Fig. This arises from the increased effect of system- 
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atic uncertainties in the event selection as the amount of 
background in the search window increases and the abil- 
ity to distinguish GRB neutrinos from background events 
becomes correspondingly more important. 
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FIG. 3. 90% CL Feldman-Cousins[13] upper limit (fluence 
normalization at 642 TeV, the first peak of the expected spec- 
trum) set by the model-independent analysis in each time 
window for an E~ 2 and for the Guetta et al. spectrum. Sys- 
tematic errors on the number of events are included. There is 
an additional t\% uncertainty on the effective area and thus 
on the right-hand axis. The three sharp peaks between 2000 
and 7000 seconds are caused by three low-energy neutrino 
events consistent with the atmospheric background. 



While the specific neutrino-flux predictions of the fire- 
ball model provided by Waxman and Bahcall|2J and by 
Guetta et al.Q are excluded (90% confidence) by this 
work, we have not yet ruled out the general picture of 
fireball phenomenology. The neutrino flux we compute 
for GRBs is determined by the flux of protons acceler- 
ated in the fireball, and by the fraction of proton energy 
transferred to charged pions The proton flux can be 
chosen either such that the energy in gammas and pro- 
tons is equal or set to the flux of cosmic rays above 10 18 
eV, with similar results. is determined largely by as- 
suming protons are accelerated, in conjunction with the 
observed low optical thickness of the source. Due to un- 
certainties in the bulk boost factor and internal structure 
of the shocks, / w may range from 10 - 30%|L3[, causing 



an uncertainty of about a factor of 2 on our calculation 
of the flux, which used w 0.2. Future observations 
by IceCube will push our sensitivity below the level of 
this theoretical uncertainty on and allow direct con- 
straints on acceleration of protons to ultra-high energies 
in Gamma Ray Bursts. 
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